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METHOD FOR B I SPECTRAL PICKING OF ANELLIPTICAL NMQ 
CORRECTION PARAMETERS 

The field of this invention is seismic prospecting. 
The invention relates particularly to processing of 
seismic traces in a common midpoint gather. 

The invention is more precisely related to a method 
5 of determining velocity V and anellipticity y\ parameters 
necessary to perform processing including NMO correction 
of seismic traces. 

Seismic prospecting usually consists of emitting 
seismic waves in the subsoil using one or several seismic 
10 sources, and making surface records of seismic data 
corresponding to seismic waves reflected on geological 
interfaces in the subsoil as a function of time (also 
called reflectors) using receivers (also called geophones 
or hydrophones) depending on whether the prospecting is 
15 being done on land or at sea) and then processing these 
data to extract useful information about the geology of 
the subsoil . 

The record of the seismic energy produced by each 
receiver during data acquisition is called a seismic 
2 0 trace. 
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One conventional technique for seismic prospecting 
is multiple coverage, in which sources and receivers are 
arranged such that several seismic traces are grouped 
together at the same midpoint (in other words the point 
5 at equal distance between the source and the receiver for 
a given trace) . 

Although seismic traces contain useful information 
about seismic reflections and geology of. the subsoil, 
they also contain noise components. 

10 One of the first purposes of the processing of 

seismic data is to eliminate, or at least attenuate, 
these unwanted noise components such that the useful 
information can be clearly identified and interpreted. 

One conventional method used to attenuate these 

15 noise components is Common MidPoint (CMP) gather. Traces 
with the same midpoint are then grouped as a function of 
the distance between the source and the receiver (called 
the offset) . 

In general, representation of seismic data in images 
2 0 requires the use of processing including: 

- a so-called TZO (Transform to Zero Offset) 
operation designed to compensate for the NMO effect of 
paths by correcting the arrival time of reflections to 
the arrival time of traces with zero offset, 
25 " and a migration operation designed to reproduce 

the correct shapes of the geological interfaces. 

Although these TZO and migration operations are 
usually carried out one after the- other, they can be 
carried out jointly. This is the case particularly when a 
30 time migration is done before stack (PSTM - Pre-Stack 
Time Migration) . 
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In a simplified manner, the TZO operation simulates 
acquisition of seismic data by sources and receivers 
placed at the common mid-point. 

The purpose is to add records illuminating the same ' 
5 point in the subsoil to increase the signal to noise 
ratio and the primary reflections to secondary 
reflections ratio, and thus benefit from the advantages 
of "multiple coverage" . 

In order to make a zero offset image, a so-called 
10 obliqueness correction method or NMO (Normal Move Out) 
correction method is used. 

If it is assumed that the subsoil is horizontally 
stratified with no lateral variation of the propagation 
velocities, it can be shown that records that have the 
15 property of illuminating the same point in the subsoil 
are records that have the same midpoint . 

However, the image of a reflection in the subsoil 
arrives at variable times depending on the offset. 
Therefore, in order to stack reflections, it is necessary 
2 0 to start by correcting the various records to bring them 
all to a common zero offset reference. 

Historically, the NMO correction is based on a 
particularly simple model of the subsoil,- a homogeneous 
model with horizontal reflectors. 
25 In this model, reflections associated with each 

subsoil reflector are theoretically aligned along 
hyperbolas also called indicators, centred vertically 
above the midpoint. 

The arrival time of a reflection is then a 
30 hyperbolic function of the source-receiver offset, the 
■shortest time being obtained at zero offset. 



WO 2005/008291 PCT/IB2004/002617 

4 

In order to . stack records in each gather, the NMO 
correction straightens the hyperbolas until they are 
theoretically horizontal. 

The NMO correction is then made based on the 
5 following hyperbolic equation for the travel time t after 
reflection, associated with a source-receiver offset x 
pair : 




in which t 0 represents the zero offset travel time and V 
10 denotes the average propagation velocity of waves in the 
subsoil . 

The simplistic model mentioned above is based 
particularly on assumptions such as low angles of 
incidence and an isotropic medium. 

15 But assumptions made for the simplistic model are 

too restrictive to describe a complex medium, and in 
particular cannot be applied to the propagation of 
seismic waves in an anisotropic medium (medium in which 
the velocity of waves can vary . depending on the 

2 0 propagation direction) . 

The use of a less simplistic model is necessary, 
particularly due to: 

- the use of long streamers for acquisition in deep 
offshore., which in particular means that traces with far 

25 offsets are acquired; 

- observation of anisotropy in clayey type 
sediments . 

The generally accepted assumption consists of 
modelling an anisotropic medium as a stack of 
30 transversely isotropic layers with a vertical axis of 
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symmetry. This is known as Vertical Transverse Isotropy 
(VTI) . 

Thus, it has been proposed to determine NMO 
corrections to be made by: 
5 - introducing vertical non- homogeneity into a 

homogenous VTI medium model, as presented in the document 
" ALKHAL I FAH T. and TSVANKIN I., 1995, Velocity analysis 
for transversely isotropic media: Geophysics, 60, 1550- 
1566"; 
10 or by 

- introducing VTI anisotropy into a model with 
stratified isotropic layers, as demonstrated in document 
"SILIQI R. and BOUSQUIE N. , 2000, Anelliptic time 
processing based on an offset hyperbola approach, 70 th 
15 Ann. Internat . Mtg. : Soc . Of Expl . Geophys . , 2245-2248". 

This second approach combining vertical non- 
homogeneity and VTI anisotropy to give a new model of the 
subsoil, appears to be the best approach in most real 
cases studied. 

The following equation for correction of the travel 
time t after reflection, using an anelliptic offset 
hyperbola, associated with a source - receiver offset x 
pair, is derived from this model : 

t(v "" ■ T^'° + £ihS Equation < la) 

where V is the velocity conventionally used in seismic 
corresponding to small offsets and t, is a parameter, 
called the anellipticity parameter. 

The document "SUAUDEAU E. and SILIQI R. , 20 01, 
Anelliptic pre-stack time migration, Annual International 
Meeting, CSEG Expanded Abstracts" also proposes to 
include the NMO correction by anelliptic offset hyperbola 
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in -the equation for the path time used during the PSTM 
migration operation. 

The PSTM migration equation is conventionally 
expressed in the form of a double anelliptic offset 
hyperbola , which is the sum of two square roots (the 
Double Square Root (DSQR) equation) . 

Taking account of anellipticity, the expression for 
this constant offset equation, becomes: 




1 + 8-n ° KU + 8T lJ (1 + 8ti)V 2 K 1 + 8r l>' (1 + 8t!)V 2 

10 Equation (lb) 

where : 

- the V and 7] parameters are those mentioned above, 
. - is the coordinates of the midpoints, 

- x - Xn represents the migration aperture, 
15 - h is the half source-receiver offset, 

- t 0 is the double time for zero aperture of the 
operator . 

Note that when the aperture x - x m of the migration 

is zero, the PSTM correction equation (lb) becomes the 
20 NMO correction equation (la) . Thus, the NMO correction 

forms a special case of the PSTM migration; the case of 

PSTM migration with zero aperture. 

Therefore, finally, in order to perform seismic data 

processing including NMO correction taking account of 
25 vertical heterogeneity and VTI type anisotropy, it is 

necessary to determine the velocity V and anellipticity r\ 

parameters . 

The said parameters V and r) can conventionally be 
estimated by carrying out two passes such that: 
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- during the first pass, the distribution of 
velocities V along the time 'axis is estimated, only using 
data with near offsets; 

- during the second pass, the anellipticity r\ is 
5 estimated along the time axis, using: 

- the distribution of velocities determined 
during the first pass, and 

- all data (including data with far offsets) . 
The document entitled: 

10 "SILIQI R. 2001, Technological leap in time 

processing focuses the data throughout anisotropic media: 
First Break,. 19, No. 11, 612-618", 

has also demonstrated that the parameters V and T) can be 
estimated in a single pass, during which bispectral 
15 analyses are carried out in order to simultaneously pick 
the V and T| parameters along the time axis, using all 
data . 

However, dense analyses of NMO correction parameters 
are preferably made when the correction equation no 
20 longer depends on the time t 0 (we then talk about static 
NMO correction) 

A static correction effectively provides a means of 
offsetting all samples each forming traces of the same 
time 5t, for a given offset. 
25 Thus, when a static correction is made, the number 

of calculations to be made can be significantly reduced 
and the trace stretching phenomenon is eliminated, which 
makes the said dense analyses viable. 

At the moment, the only way to make velocity and 
30 anellipticity analyses is to use two passes, in 
particular making use of parabolic approximations of NMO 
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residuals, in order to obtain a dense picking of the V 
and r| parameters. 

In this context: 

- the first step is to estimate residual velocities, 
5 using near offset data, according to a first estimate of 

velocities; 

- the next step is to estimate anellipticity on all 
data, using velocity updates done previously. 

A mute law must also be defined to estimate residual 
10 velocities, so as to only keep data that are considered 
as being near offset data. 

However, the efficiency of the two-pass ' analysis is 
particularly sensitive to the choice of such a mute law. 

On the other hand, the anellipticity r\ is estimated 
15 mainly making use of far offset data. 

But, the anellipticity estimate made for far offset 
data is not very precise, such that the correction made 
is finally inaccurate. 

One purpose of the invention is to overcome these 

2 0 limitations and disadvantages, by proposing a processing 

including a static NMO correction which is more efficient 
and more precise. 

More precisely, the purpose of the invention is a 
dense determination of velocity V and anellipticity r\ 
25 parameters in a single pass using all available data, in 
other words making use of all offset ranges. 

To achieve * this, the invention proposes a method of 
determining the velocity V and anellipticity r\ parameters 
for processing seismic traces from a common midpoint 

3 0 (CMP) gather including an anelliptic NMO correction, 

comprising: 



WO 2005/008291 



9 



PCT/EB2004/002617 



- a preliminary step to define a plurality of nodes 
(dtn, x 0 ) , the said nodes being indicative of parameters 
dtn and t 0 representing the NMO correction for the 
maximum offset and the zero offset travel time in 

5 hyperbolic coordinates, the said preliminary step being 
followed by 

- for each node (dtn, x 0 ) defined in the preliminary 
step, steps for the NMO correction of traces in the CMP 
gather as a function of the values of the said parameters 

10 dtn and t 0 at the node considered, and for calculating 
the semblance function associated with the said NMO 
correction for the node considered; and 

- for each picked time t 0/ a step including 
determination of the maximum semblance node (dtn(t 0 ), 

15 To (t 0 )), 

- and a final step to convert the dtn(t 0 ) and x 0 (t 0 ) 
parameters so as to obtain the velocity V(t 0 ) and 
anellepticity r\ (t 0 ) laws. 

According to a first embodiment of the invention, 
20 the processing done is a' static NMO correction of seismic 
traces . 

According to a second embodiment of the invention, 
the processing done is a PSTM migration of seismic 
traces, the said PSTM migration including a static PSTM 
25 NMO correction of the said seismic traces. 

A preferred but non-limitative aspect of the process 
according to the invention relates to the definition of 
the dtn and x 0 parameters relative to the velocity V and 
the anellipticity T| so as to make static NMO corrections, 
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according to dtn = _2ll_t + I _Ifi_ + _JSns* and _ = 



tr 



1+8ti 

Other aspects, purposes and advantages of this 
invention will become clearer after reading the following 
detailed description with reference to the attached 
Figures in which: 

- Figure la shows the anelliptic offset hyperbola 
used to make the NMO correction and illustrates the 
meaning of the parameters t 0 and dtn; 

- Figure lb represents the DSQR equation of the 
double anelliptic offset hyperbola of the PSTM migration 
and illustrates the meaning of the t 0 and dtn parameters; 

- Figure 2 illustrates the effect of the t 0 
parameter of the reflection curve corrected by the NMO 

15 correction; 

- Figures 3a and 3b represent the analysis volume 
(t 0/ dtn, t 0 ) in which the bispectral picking of the dtn 

•and T 0 parameters according to the invention is made; 

- Figure 4 combines each of (V, V m ) and (dtn, t 0 ) 
20 approaches, representing their corresponding bispectral 

analysis panel; 

- Figure 5 shows the correspondence between the 
(dtn, T 0 ) and (V, V^) pairs of NMO correction parameters; 

- Figure 6 shows a CMP gather of true seismic traces 
2 5 before the NMO correction and the bispectral picking of 

the (dtn, To) parameters corresponding to this gather of 
traces ; 
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- Figure 7 shows the semblance function and velocity 
and anellipticity functions deduced from the picking of 
the dtn and t 0 parameters in Figure 6. 

- Figure 8a is a flowchart representing the steps of 
5 a first particular embodiment of the invention, namely 

determination of the V and r| parameters to make an NMO 
correction; 

- Figure 8b is a flowchart representing the steps of 
a second particular embodiment of the invention, namely 

10 determination of V and r| parameters to do a PSTM 
'migration; 

- Figure 9 is a diagram illustrating the different 
operations carried out in order to determine the V and rj 
parameters to make a PSTM migration. 

15 In general, the method according to the invention is 

a method for processing records of variable offset 
seismic traces, and this processing uses the recorded 
seismic traces to build common midpoint (CMP) trace 
gathers, and the traces in each gather are subjected to 

20 an NMO correction. 

In particular, the method according to the invention 
determines the velocity V and anellipticity r| parameters 
to perform processing including such an NMO correction of 
the seismic traces of a CMP gather. 

25 Por example, the said processing may be: 

- an anelliptic NMO correction of the seismic 
traces ; 

- an anelliptic PSTM migration which, as we have 
seen above, jointly applies TZO and migration operations 

3 0 (in the following we will refer to a PSTM NMO 
correction) . 
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The following description is more specifically 
applicable to the NMO correction. However, considering in 
particular the two particular embodiments of the 
invention that will be described below, we can see that 
this description is equally applicable to any processing 
including an NMO correction, particularly processing 
including a PSTM NMO correction. 

Two new parameters dtn and t 0 are considered in 
order to specifically determine V and tj: 

- To that represents the zero offset travel time in 
"hyperbolic coordinates" (see Figure 1) , 

T ° = l + &n 11 Equation (2) 

and, 

- dtn that represents the NMO correction for the 
15 largest offset x ma x (see Figure 1), 

dtn = t x=Xmax - t Xe0 namely: 



10 



dtn = -23-to+. fl 



l o 



x ma. Equation (3) 



(1 + 8ti)V 2 

. ■ It is important to note that dtn is defined relative 

to the velocity V and the anellipticity t\, while x 0 is a 
20 perfectly anelliptical parameter defined relative to the 
anellipticity r\, independent of V. 

Accordingly, the velocity V and anellipticity r| 
parameters may be calculated in accordance with equations 
(2) and (3) , using the following conversion equations (4) 
25 and (5) : 



— Equation (4) 



dtn(dtn + 2T 0 )^- 



o 
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and ri = - 
1 8 



— -1| Equation (5) 

\ T * J 



Figure la shows the anelliptical offset hyperbola 
5 used to make the NMO correction and illustrates the 
meaning of the x 0 and dtn parameters . 

The said "hyperbolic coordinates" are shown in this 
Figure 1. Their origin on the time axis is taken at the 
intersection of the said time axis with the asymptote 
10 tangent to the said far offset offset hyperbola (see 
equation (la) ) . 

Using the (dtn, t 0 ) parameters, equation (la) for 
the anelliptical offset hyperbola becomes: 



Hi dtn(dtn + 2x 0 ) 2 „ ... ,^ v 
t = t 0 - T 0 + TJ5+ - — ^x' Equation (6a) 

V X max 

15 The (dtn, t 0 ) parameters defined for the velocity V 

and anellipticity tj are thus used ' to make the NMO 
correction CORRnmo = t - t 0 to be applied to offset x 
traces independent of t 0 : 



\o dtn(dtn + 2T n ) -> 
CORRnmo(x) = -To + |To+ 2 — x Equation (7a) 

V X max 

20 Therefore, this is a static NMO correction. In other 

words, the data recorded on a given offset trace will all 
be corrected in the same way for a (dtn, x 0 ) pair, 
independently of the time at which these data were 
acquired. 

25 Consequently, the estimate of the velocity and 

anellipticity parameters is not disturbed by stretching 
of traces usually observed when dynamic NMO corrections 
are made . 
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Similarly using the two parameters (dtn, t 0 ) , 
equation (lb) for the DSQR double offset hyperbola for 
the anelliptic PSTM migration becomes: 



t - t 0 - t 0+ p , dtn ( dtn + 2T oX x -^+h) 2 + 

V ^ X max 



^ dtn(dtn + 2x 0 )(x - x m - h) s 
— + %± 52 Equation (6b) 

x max 



Figure lb shows the DSQR double offset hyperbola of 
the PSTM migration defined with the (dtn, r 0 ) parameters. 

In the context of the PSTM migration, x^ represents 
the maximum offset and aperture of the migration. 

The (dtn, t 0 ) parameters defined relative to the 
velocity V and anellipticity T) are thus used to make the 
PSTM NMO correction CORR PSTM = t - t 0 to be applied to the 
offset x traces independent of t 0 (equation 7b) ) : 

Corr PSTM (x) = -To + R , dtn(din + 2r 9 )(x-x + /^ + /r„ 2 t tftnfrfoi + 2t 0 )(x - X - h) 2 

V. 4 V4 + 

Therefore, this is a static PSTM correction. In 

other words, the se of aperture samples x - x* of an 

-iso-offset" cube are offset by the same time for a given 
(dtn, t 0 ) pair. 

By analysing the (dtn, t 0 ) parameters in several 
picking times, the method according to the invention can 
be used in particular to determine the (V, r\) parameters 
necessary for processing including an anelliptic NMO 
correction of the traces of a CMP gather. 

The said process includes the steps presented be low- 
in a simplified manner. 
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During, the preliminary step, an analysis volume is 
defined including several nodes (dtn, t 0 ) • 

The following are performed for all nodes in this 
volume : 

5 - firstly, according to equation (7a) , the static 

NMO correction of traces in the studied CMP gather, as a 
function of the values of the dtn, t 0 parameters at the 
node considered, the said static correction being valid 
for any picking time; 

10 - secondly, the semblance as a function of time 

associated with the correction made in the previous step 
is calculated. 

Finally, the following are determined for each time 
in the several picking times: 

15 - firstly, the node (dtn, x 0 ) used to make an 

optimum correction, for example with regard to the 
semblance criterion (this type of criterion typically 
* being used in seismic processing to "measure the 
horizontality" of reflection curves and to determine the 

20 reliability of the picking) ; 

- secondly, the values" of the dtn, t 0 parameters at 

the said maximum semblance node are converted to values 

of the velocity V and anellipticity r\ parameters at the 

said picking time considered. 
25 The velocity law (in other words all (picking time, 

V)) pairs and the anellipticity law (all (picking time, 

r\) ) pairs are thus defined. 

Finally, the NMO correction for all seismic traces 

may be made using, these velocity V and anellipticity r\ 
30 laws in equation (la) for the anelliptic offset 

hyperbola. 
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Similarly, and as will be described in detail later, 
the (dtn, t 0 ) parameters of the static PSTM migration 
(see equation 7b) ) may also be checked. The velocity V 
and anellipticity r\ laws are then determined and may be 
5 used to make the PSTM migration in the DSQR equation (lb) 
for the double anelliptic offset hyperbola. 

Figure 2 illustrates the effect of the parameter x 0 
(and therefore the anellipticity r\ according to equation 
(2) on the curvature residues after the NMO correction. 
10 Note that the vertical time scale of the curve in 

Figure 2 is exaggerated so that this effect can be 
clearly understood. 

Three curves are shown in Figure 2 , for which the 
dtn parameter is fixed to the correct value and the 
15 parameter t 0 is assigned different values. 

The central curve shows the case in which t 0 is 
equal to its correct value t 0 i, in other words when the 
corresponding anellipticity t) x is equal to its true value 
iltrue. As expected, the corrected reflection curve is then 
2 0 horizontal. 

The upper curve shows the case in which t 0 is equal 
to a value x Q2 less than its correct value t 0 i, the 
corresponding anellipticity r\ 2 being greater than its 
true value Titme- 
25 The lower curve represents the case in which x 0 is 

equal to a value t 03 greater than its true value x 01 , the 
anellipticity r| 3 being less than its true value Ti true . 

Note that these lower and upper curves show that the 
"horizontality" of the corrected reflection curve is 
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acceptable at near offset (x « 0) and at far offset (x « 

On the other hand, significant curvature residues 
are observed when the offset x no longer tends towards 
5 one of these limiting values 0 and Xma*. In particular, 
particularly significant residues are observed for an 
offset x centred in the middle of the offsets range. 

For example, and as is shown in Figure 2, when x 0 is 
equal to x Q2r an RMOx 0 2 residual correction must be made 
10 to the offset traces XrmoT 0 2. Similarly, when x 0 is equal 
to x 03/ a residual correction RMOt 0 3 must be made to 
offset traces Xrm O ^03 • 

Due to these significant curvature residues, the 
offsets range can be used almost in its entirety in order 
5 to determine the anellipticity t\. 

Therefore, setting parameters in (dtn, x 0 ) for the 
NMO correction make it possible to use available data for 
all offsets (x between 0 and Xwax) in determining the 
anellipticity r\ . 
0 As has already been mentioned, thi£ is not. the case 

for NMO corrections for which parameters have been set 
with V and r\ for which the estimate of the anellipticity 
r| is made essentially using far offset data. 

Thus, the effect of the new anelliptic parameter t 0 
5 is distributed on all offsets, unlike the anellipticity r\ 
that only affects far offsets. Therefore, the "behaviour" 
of the parameter t 0 gives a better constraint on 
anellipticity values. 
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As already mentioned above, the optimum dtn and t 0 
parameters are determined within a 3D analysis volume 
(to# dtn, x 0 ) . 

Several nodes (dtn, x 0 ) , in other words several 
5 pairs of dtn, x 0 parameters for which values are known, 
are considered in the said analysis volume. 

Nodes are usually regularly spaced from each other, 
by an increment Adtn on the dtn axis and an increment At 0 
on the t 0 axis, 

1° Minimum values dtrw*, x 0m i n , t 0 m±n and maximum values 

dtn maX / x 0 max# tomax of the x 0 and to respectively provide a 
means of defining the limits of the said analysis volume. 

Advantageously, plausible values of the velocity 
parameter V and the anellipticity parameter r\ may be used 

15 to define a corridor [dtn min (t 0 ) , dtn max (t 0 ) ] , [x 0m in(t 0 ) , 
toroax(to)] inside the said analysis volume. 

This corridor restricts the analysis volume and 
therefore the number of nodes (dtn, x 0 ) that have to be 
considered to determine the optimum pair, (dtn, x 0 ) . 

2 0 Although use of the said corridor is beneficial for 

the efficiency of the process according to the invention, 
it also provides a means of constraining solution towards 
the right phenomena, without needing to consider 
incompatible (dtn, x 0 ) (and therefore V, ti) ) pairs, for 
25 example related to multiple reflections or miscellaneous 
interference phenomena. 

Figures 3a and 3b represent the said analysis volume 
(t 0/ dtn, T 0 ) in the context of an example of true seismic 
data processing during an NMO correction done according 

3 0 to the process according to the invention. 
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Figure 3a represents three 2D panels a, b and c of 
the analysis volume: 

- panel a is a panel (dtn, x 0 ) with constant t 0 ; 

- panel b is a panel (dtn, t 0 ) with constant x 0 ; 
5 - panel c is a panel (t 0 , t 0 ) with constant dtn. 

Figure 3b diagrammatically shows the 3D analysis 
volume (t 0 , dtn, t 0 ) and three intersections with this 
volume along three planes with constant t 0/ t 0/ dtn 
respectively, each of these intersections being projected 

10 onto the corresponding side a, b or c in Figure 3a. 

Panels (dtn, x 0 ) with constant t 0 (panel a in the 
above example) are panels in which the bispectral picking 
of the dtn and x 0 parameters is made, for example 
according to the maximum NMO correction semblance 

15 criterion, for the considered picking time t 0 . 

The [dtn min (t 0 ), dtn ma x(t 0 >] [T 0m in (t 0 ) Tomax (t 0 )] 

corridor mentioned above for the effective analysis 
inside the analysis volume is " also shown in Figure 3a 
(lighter zone on each panel) . 
20 The seismic data sampling time step At 0 defines the 

difference between two successive picking times (times 
for which the maximum associated semblance node (dtn, t 0 ) 
is determined) and therefore the number of bispectral 
picking panels to be considered. 

2 5 Automatic picking also provides a means of 

extracting the dtn, r 0 parameters at a density that is 

greater when the At 0 increment between picking times is 
small . 

Sampling of the analysis parameters dtn and t 0 is 

3 0 directly related to the resolution of the seismic 
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exploration; dtn and x 0 effectively have the same 

dimensions as seismic records. 

The systematic search for the maximum semblance, 

conventionally known in itself, provides a means of 
5 determining the (dtn, x 0 ) pair providing best focus, for 

a given picking time t 0 . 

Parabolic interpolations about the values of nodes 

(dtn, x 0 ) can also provide a means of evaluating values 

of the dtn, x 0 parameters between the different nodes 
10 that have actually been picked. And in particular this 

type of evaluation makes determination even more precise 

(in contrast to determination limited to nodes in the 

corridor) of the dtn, t 0 parameter pair maximising the 

semblance function. 
15 Finally, the velocity V and anellipticity r\ 

parameters are determined, always for the picking time t 0 

considered, by using the above mentioned conversion 

equations (4) and (5) . 

Figure 4 shows the comparison between the two 
20 approaches (V, Van) and (dtn, . x 0 ) / by representing the 

bispectral picking panel for each, for a given picking 

time . 

The figure at the right illustrates the conventional 
approach (V, V an ) for which the two axes are the velocity 
25 axes (the anellipticity t] being related to the ratio of 



these two velocities according to r\ = ^ 



v^ 4 J 



The left figure illustrates the approach (dtn, x 0 ) 
according to the invention for which the two axes are 
time axes. 
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It is important to note from the study in Figure 4 
that the dtn and t 0 parameters appear to be decorrelated. 
This " decor relation" is striking when the two approaches 
are compared, spreading of the (dtn, x 0 ) spectrum 
5 actually being much narrower than spreading of the (V, 
Van) spectrum. 

Therefore the picking made in the context of the 
(dtn, t 0 ) approach according to the invention is more 
precise than the approach carried out in conventional 
10 methods . 

Furthermore, this decorrelation enables filtering of 
the dtn and t 0 pickings separately, while keeping NMO 
corrections. This is not the case for the V and T| 
parameters for which a reduction in one of the parameters 
15 must necessarily be compensated by an increase in the 
other and vice versa. 

And due to interpolations and individual filtering 
of parameters according to the invention, dtn and x 0 , it 
is then possible to perform simultaneous interpolation 
20 and filtering of the standard NMO correction parameters V 
and *n . 

Figure 5 shows non- linear correspondence according 
to equations (4) and (5) mentioned above, between the 
pair of time parameters (dtn, t 0 ) and the pair of 
25 velocity parameters (V, V^) . 

Figure 6 shows, from left to right: 

- a CMP gather of real seismic traces before NMO 
correction; 

- the picking of the dtn parameter corresponding to 
3 0 this gather of traces; 
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w the picking of the x 0 parameter corresponding to 
this gather of traces. 

The straight line t 0 = t 0 on the picking of the 
parameter x 0 at the right in Figure 6, corresponds to the 
5 purely hyperbolic reflection curves. 

From right to left, Figure 7 represents the 
semblance function and the velocity V and the 
anellipticity r\ functions (see equation (4) and (5) ) 
deduced from the picking of the dtn and t 0 parameters 
10 shown in Figure 6. 

On this example of real seismic data processing, it 
can be seen that the values of V and ri obtained generally 
corresponding to a semblance of more than 40%. 

The following description contains details of two 
15 particular embodiments of the invention. 

The first of these modes relates to a process for 
determining optimum parameters to make an anelliptic NMO 
correction to the traces of a CMP gather (see the 
different steps shown in the flowchart in Figure 8a) . 
20 With reference to Figure 8a, this first embodiment 

• includes an initialisation- step la during which the 
following operations are carried out in sequence: 

- determination of the limits of the analysis volume 

[dtrimin/ dtn max ] t^Omin/ ^Omax] [tomin/ tomax] ; 

25 - calculation of NMO corrections CORRnmo (equation 

(7a)) for all offsets and for all nodes (dtn, x 0 ) 
included in the analysis volume; 

- delimitation inside the analysis volume of the 
corridor [dtn rain (t 0 > , dtxw(to)] [xomin (to) ^omax <t 0 ) ] of 

30 plausible velocity and anellipticity values. 
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Once the initialisation step la has been carried 
out, a step 2a to calculate velocity V(t 0 ) and 
anellipticity r) (t 0 ) laws is performed for each gather of 
CMP traces . 
5 This step 2a comprises: 

- a first operation 3a carried out for each node 
(dtn, x 0 ) in the corridor defined in the initialisation 
step la, during which the following operations are 
carried out in sequence for each picking time t 0 : 

1° - application of static NMO corrections CORRnmo 

for all offsets along the corridor, precalculated during 
the initialisation step la; 

- calculation of the semblance function on data 
corrected along the corridor using a time window 

15 appropriate to the dominant wavelet; 

- summation ("stack" . calculation) of data 
corrected along the corridor (only near offset data can 
be used advantageously for this purpose) ; 

- a second operation 4a carried out for each picking 

2 0 time t 0 (the said times being at intervals of At 0 between 

rtomin/ tomax] during which the following operations are 
carried out : 

- search for maximum semblance in the [dtn^inCto) , 
dtn^axfto)] [xomin (t 0 ) T 0m ax (t 0 )] corridor of the bispectral 

2 5 panel (dtn, t 0 ) ; 

- check the fact that the position in (dtn, x 0 ) at 
the maximum semblance corresponds to a summation extreme 
value for the same values dtn and x 0 ; 

- creation of the dtn(t 0 ), T 0 (t 0 ) and semblance 
30 ' (t 0 ) series; 
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- a third operation 5a designed to select and adjust 
the pickings obtained, during which the following 
operations are carried out: 

- increasing sort of the semblance series (t 0 ) / 

5 - validation of pickings dtn and x 0 for which the 

distance in time to the highest semblance pickings is 
greater than a predefined value; 

- adjustment of picking dtn and t 0 values 
validated by parabolic interpolations using surrounding 

10 values; 

- retention of picked, validated and adjusted 
values if it is possible to calculate the Dix interval 

• velocities with the pickings with the highest 
similarities . 

15 - a fourth operation 6a designed to use equations 

(3) and (4) to convert picked, validated and adjusted and 

retained values of dtn and t 0 during operation 5a, into 

velocity V and anellipticity r\ laws. 

The velocity V and anellipticity ti laws as a 
2 0 function of time are thus perfectly determined. The 

anelliptical NMO correction of seismic traces in the CMP 

gather can thus be done precisely. 

The second particular embodiment of the invention 

relates to a process for determination of optimum 
25 parameters for anelliptical PSTM migration of traces in a 

CMP gather. 

This second embodiment may be included in a 
generalisation of the first mode discussed above. 

As has already been demonstrated, the use of the 
30 (dtn, x 0 ) parameters, enables static PSTM NMO corrections 
(see equation (7b)>. 
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In the context of a PSTM NMO correction, this 
application has the same advantages as previously 
discussed for the NMO obliqueness correction. 

More precisely, it will be noted that the first 
5 embodiment is only a special case of the second 
embodiment corresponding to the case of a zero migration 
aperture . 

With reference to Figure 8b, the second embodiment 
includes an initialisation step lb during which the 
10 following operations are carried out in sequence: 

- determine the limits of the analysis volume 
[dtn min , dt nmax ] , [t 0mia( Tomax ] , [tomin/ tomax] . 

- calculate NMO corrections CORR PSTM (equation (7b)) 
for all nodes (dtn, t 0 ) included in the analysis volume 
and for all migration offsets inside the migration 
aperture ; 

- delimitation of plausible velocity and 
anellipticity values inside the analysis volume of the 
corridor [dtrWto), dtn max (t 0 )3 [T 0min (t 0 ) T 0max (t 0 ) ] . 

Once this initialisation step lb has been done, the 
said first embodiment creates' a step 2b to calculate the 
velocity V(t 0 ) and anellipticity * (t 0 > laws for each 
gather of CMP traces. 

This step 2b comprises: 

- a first operation 3b done for each node (dtn, x 0 ) 
of the corridor defined during the initialisation step 
la, during which the following are carried out in 
sequence : 

- for each offset class, operations for: 

- application on all midpoints inside the 
migration aperture along the corridor, of static 
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corrections C0RR PS t M precalculated during the 
initialisation step lb; 

- summation of corrected midpoints along the 

corridor . 

- for each picking time t 0 , operations to: 

- calculate the semblance function on 
corrected data along the corridor using a time window 
appropriate to the dominant wavelet; 

- summation (stack calculation) of data 
corrected along the corridor (only data with small 
offsets can be used advantageously for this purpose) : 

- a second operation 4b carried out for each picking 
time t 0/ similar to the operation 4a described above, to 
create the dtn(t 0 ), x 0 (t 0 ) and semblance t 0 series; 

15 ~ a third operation similar to operation 5b 

described above, designed to select and adjust the 
pickings obtained; 

- a fourth operation 6b similar to operation 5a 
described above, designed to convert values of dtn, t 0 

20 into velocity V and anellipticity r\ laws. 

The velocity V and anellipticity t\ laws as a 

function of time are thus perfectly determined. And the 

PSTM migration of seismic traces in the CMP gather may 

thus be done accurately. 
25 Figure 9 illustrates the second embodiment of the 

invention that has just been described. 

Seismic data are initially grouped into iso-offset 
cubes . 

As has just been described, NMO corrections CORR PSTM 
are applied (see first operation in step 3b) for each 
offset class (namely for each iso-offset cube) on all 
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midpoints. The arrow marked with the MIG label in Figure 
9 illustrates this, operation. 

The midpoints thus corrected are then summated 
during a second operation in step 3b, the arrow labelled 
5 STACK in Figure 9 illustrating this operation. 

These two MIG and STACK operations are specific to 
the second embodiment of the invention (PSTM migration) . 

The following operations are performed in the 
context of each of the two embodiments discussed (the NMO 
10 obliqueness correction, as we have described, 
corresponding to the special case of a PSTM migration 
with zero migration aperture) . 

The semblance is then calculated for each picking 
. time (second operation in step 3a for the NMO correction 
15 only, third operation in step 3b for the PSTM migration) , 
. the arrow labelled 11 semblance" in Figure 9 illustrating 
this operation. 

Obviously, the MIG, STACK and semblance calculation 
operations are implemented for each node (dtn, x 0 ) . 
20 The "automatic bispectral picking" mentioned in 

Figure 9 corresponds to picking of maximum semblance 
parameters (dtn, x 0 ) for each picking time t 0 (operations 
4a, 4b in Figures 8a, 8b respectively) . The arrow 
labelled MAX illustrates the search for the maximum 
25 semblance. 

Finally, the pickings of the picked parameters (dtn, 
t 0 ) are converted into velocity V and anellipticity y\ 
laws (operations 6a, 6b in Figures 8a, 8b and arrow 
labelled CONV in Figure 9) . 



